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The complete knowledge of the nucleon spin structure at leading twist requires also addressing the 
transverse spin distribution of quarks, or transversity, which is yet unexplored because of its chiral- 
odd nature. Transversity can be best extracted from single-spin asymmetries in fully polarized Drell- 
Yan processes with antiprotons, where valence contributions are involved anyway. Alternatively, 
in single-polarized Drell-Yan the transversity happens convoluted with another chiral-odd function, 
which is likely to be responsible for the well known (and yet unexplained) violation of the Lam- Tung 
sum rule in the corresponding unpolarized cross section. We present Monte-Carlo simulations for 
the unpolarized and single-polarized Drell-Yan pp'^' — > ii~ X at different center-of-mass energies 
in both configurations where the antiproton beam hits a fixed proton target or it collides on another 
proton beam. The goal is to estimate the minimum number of events needed to extract the above 
chiral-odd distributions from future measurements at the HESR ring at GSI. It is important to study 
the feasibility of such experiments at HESR in order to demonstrate that interesting spin physics 
can be explored already using unpolarized antiprotons. 

PACS numbers; 13.75.-n,13.85.-t,13.85.Qk,13.88.-|-e 



I. INTRODUCTION 

When studying the spin structure of hadrons, the best tool is represented by the spin asymmetry, defined by the 
ratio between the difference and the sum of differential cross sections obtained by flipping the spin of one of the 
involved polarized hadrons. For hadronic beams, spin asymmetries are known since a long time; historically, the first 
one was obtained at FERMILAB, where an anomalous large transverse polarization of the A produced in proton- 
nucleon annihilations was measured [J , surviving even at large values of the A transverse momentum. More recently, 
similar anomalously large asymmetries have been observed by the FNAL E704 ^ and STAR collaborations in 
inclusive pion production from collisions of transversely polarized protons. Such asymmetries require a nonvanishing 
imaginary part in the off-diagonal block of the fragmentation matrix of quarks into the detected hadron, which is 
forbidden in perturbative QCD at leading twist '?|. Soon after the first FERMILAB observation, a pioneering work 
appeared 5] about the possibility of having leading-twist asymmetries in fully polarized Drell-Yan processes, based 
on a new partonic function describing the quark transverse spin distribution in a transversely polarized hadron: the 
transversity. However, this work was basically ignored (and the observed asymmetries left unexplained) upon the 
prejudice that transverse-spin effects have to be suppressed. Only almost a decade later, it was realized that the 
transversity is connected to a helicity fiipping mechanism but it is still a leading-twist distribution, necessary to 
complete the knowledge of the nucleon spin structure 0,0 (for a review, see also Refs. 

As such, the transversity is not diagonal in the parton helicity basis; since at leading twist helicity and chirality 
coincide, it is usually mentioned as a chiral-odd function. But in the transverse spin basis it is diagonal and can 
be given a probabilistic interpretation. To be precise, in a transversely polarized hadron it describes the probability 
difference for having a transversely polarized parton with spin parallel or antiparallel to the one of the parent hadron. 
The chiral-odd feature brings in negative and positive consequences. The transversity needs a chiral-odd partner 
in order to be measured in a (chiral-even) cross section. At leading twist, the partner can be the transversity 
itself (for the corresponding antiquark in the aforementioned polarized Drell-Yan) or another chiral-odd unknown 
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fragmentation function in a semi-inclusive process. In particular, the transversity is not accessible in inclusive Deep- 
Inelastic Scattering (DIS) Therefore, in the Quark Parton Model (QPM) description of DIS, it does not have a 
counterpart at the level of structure function. Even if it depends on spin, it is not related to some partonic fraction of 
the hadron spin; the related twist-2 tensor operator a^'^j^ is not part of the hadron full angular momentum tensor, 
and its first moment, the tensor charge, has a nonvanishing anomalous magnetic moment: contrary to the nonsinglet 
axial charge, from a starting input scale the tensor charge unavoidably evolves to zero j.8]. 

If the hadron target has spin ^ , there is no transversity for gluons because of the mismatch in the change of helicity 
units. Hence, during evolution the quark transversity decouples from any radiative gluon. Moreover, it changes sign 
under charge conjugation, which implies that again during evolution it does not receive any contribution from charge- 
even structures like the qq pairs from the Dirac sea. In summary, the very same features that make the transversity 
quite an elusive object, are responsible for its very peculiar behaviour. Its evolution is similar to the one of a nonsinglet 
structure function; it describes the spin distribution of a valence quark while "switching off the contribution of gluons 
and of the vacuum" j^j. 

The growing interest in the transversity reflected in a rich experimental program, that led to some first recent 
measurements with hadronic TT] and leptonic beams [l^ lisl IT3 ITsf. However, the field is confronted with rapid 
developments that should bring new upgraded and precise results in the near future |l6j . Transverse spin is intimately 
connected with an explicit dependence of the parton distribution on its intrinsic transverse momentum. Including this 
explicit dependence, the list of parton distribution (and fragmentation) functions is very rich already at leading twist, 
containing also other chiral-odd objects 17]. While on the theoretical side the discussion is focussed on extending 
the universality proof to such function s [ip. 11^ pdi llll |2^ . new azimuthal asymmetries are being deviced and 
measured both on the phenomenological |23ll24ll25| | and experimental [2^ side, that can shed light on new underlying 
mechanisms, ultimately related to the orbital motion of the quarks in the parent hadron. 

In principle, the aforementioned polarized Drell-Yan is the most convenient process to measure the transversity at 
leading twist, because it does not involve other unknown functions. However, a Next-to-Leading (NLO) simulation 
has shown that the QCD evolution and the Soffer inequality (if forced to be fulfilled at any scale during evolution) 
make the resulting spin asymmetry very small '27', '2^. Moreover, the latter involves the transverse spin distribution 
of an antiquark in a transversely polarized proton, which is not expected to be large. This difficulty can be overcome 
at the foreseen HESR ring at GSI, where (polarized) antiprotons will be produced [2^ Isollsil ls^ Iss f : in the process 
pTpT _> l+l~X, only valence (anti)quark distributions are involved without any further suppression. 

But interesting information on the partonic spin structure of hadrons can also be inferred by using unpolarized 
antiproton beams only. In fact, in single-polarized Drell-Yan like pp^ l^l^X the polarized part of the cross section 
contains at leading twist three terms of which two produce interesting azimuthal spin asvmmetri es l34l . The first one 
involves the convolution of the standard unpolarized distribution /i with the Sivers function |35j , which describes 
how the distribution of unpolarized quarks is affected by the transverse polarization of the parent proton. The very 
same /f^ appears at leading twist in semi-inclusive DIS processes like Ip^ I'ttX Q or pp^ — > ttX [s^I, and it is 
responsible of an azimuthal asymmetric distribution of the detected pions depending on the direction of the target 
polarization, the socalled Sivers effect. A measurement of a nonvanishing asymmetry would be a direct evidence of 
the orbital angular momentum of quarks [16j . 

The second term in the single-polarized Drell-Yan cross section has a different azimuthal dependence and it involves 
the transversity hi convoluted with the chiral-odd distribution hj^ 34], which describes the influence of the quark 
transverse polarization on its momentum distribution inside an unpolarized parent hadron. Extraction of the latter 
is of gre at irnportance, because hj^ is believed to be responsible for the well known violation of the Lam- Tung sum 
rule |38l l39l lloj , an anomalous big azimuthal asymmetry of the corresponding unpolarized Drell-Yan cross section, 
that neither NLO QCD calculations |^, nor higher twists or factorization-breaking terms in NLO QCD 42, -13. lil] 
are able to justify. 

In this paper, we will present numerical simulations for the unpolarized and single-polarized Drell-Yan process 
pp(T) _> n+jj^X at several center-of-mass energies using both configurations where the antiproton beam hits a fixed 
proton target or it collides on another proton beam. The goal is to estimate the minimum number of events needed to 
make the extraction of hi and hi feasible in future measurements at the HESR ring at GSI. In Sec.m the kinematics 
and the general cross section formulae are reviewed. In Sec. IIIII details about the numerical simulation are given and 
results are presented in Sec. IIVI Finally, some remarks and conclusions are drawn in Sec. 

II. THE GENERAL FRAMEWORK 

In a Drell-Yan process, two (polarized) hadrons Hi and H2, carrying momentum Pi and P2 with Pi(2) — -^^1(2) 
(and spin Si, S2, such that = ^1(2) ' '^'1(2) — 0), annihilate into a lepton I and an antilepton / with momenta 
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FIG. 1: The leading-twist contribution to the Drell-Yan process. 



ki and k2 (with k^^^) ^ 0)- The available initial squared energy in the center-of-mass (cm) frame is s = (Pi + P2) ■ 
The time-like momentum transfer defines the hard scale and it is directly related to the invariant mass of the final 
state, i.e. = = {ki + ^2)^ > 0. Therefore, it is also indicated as = AP; we will use both notations, in the 
following. The DIS regime is defined by the limit Q^, s 00, but keeping the ratio < r = Q^/s < 1 limited. In 
this regime, the Drell-Yan process can be factorized into an elementary annihilation process, convoluted with two soft 
functions describing the distribution of the two annihilating partons inside the corresponding hadrons. 

If M is constrained inside a range where the elementary annihilation can be safely assumed to proceed through 
a virtual photon converting into the final 1 1 pair, while exploring the largest available range for r, the dominant 
contribution at leading order is depicted in Fig. where the blobs represent the correlation functions for the 

annihilating antiparton (labelled "1") and parton (labelled "2"), respectively: 

^P2;P2,S2) = J ^e^P^-^P2,S2\m^{^)\P2,S2) . (1) 

When — > oo, the hard scale selects a light-cone dominant component for the hadron momenta, namely and 
P2 . Then, the parton momenta Pi{2) are approximately aligned with the corresponding Pi{2) and are given by the 
light-cone fractions 



P+ 2Pi-g'^' P- 2P2-q 



^1 = — - ^715— ' - — ^ Tnr-7. ' ^ ^1(2) ^ 1 ' (2) 



while g"*" = p'^ and q~ = P2 , by momentum conservation [33 | . 

As already anticipated in Sec. ^ it is crucial to keep memory of the intrinsic transverse-momentum dependence 
inside the correlations functions of Eq. (Q. Hence, the cross section will be kept differential in q^, which is bounded 
to = Ply -|- P2y by momentum conservation; p^^ is the transverse component of the parton momentum pi with 
respect to the axis defined by the corresponding hadron 3- momentum P.^, with i — 1,2. In this context, if q^, 7^ the 
annihilation direction is not known. Hence, it is convenient to select the socalled CoUins-Soper frame [43 (see fig.|2l, 
where 



i = ^ 



Q 

XiPi X2P2 



Q Q 



(3) 



The azimuthal angles lie in the plane perpendicular to t and z. In particular, 0, , are the angles of h, S.^^ , {i — 1,2) 
with respect to the lepton plane, respectively, and fc'^ = {g'^'^ — i^f^ + z'^z'^)kii, the perpendicular component of fcf , 
with i = 1,2. Moreover, the cross section will be kept differential also in the solid angle for the lepton production. 
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i.e. dfl = 2dyd4>^ where the invariant y — /q reduces to (1 + cos0)/2 in the CoUins-Soper frame, with 9 defined 
in Fig. 121 

The hadronic tensor is obtained by taking the trace of the above correlation functions, after integrating upon the 
corresponding suppressed Hght-cone directions: 



p+=2;iPj+,P2 =X2P2 



(4) 



Using the leading-twist parametrization for $ and <i> with explicit dependence on pi^ and p2j. [4a|, respectively, the 
fully differential cross section for the unpolarized Drell-Yan process H2 ^ I IX can be written as 



da" 



dVtdxidx2dci^ 3Q 



= ^E4^(i/)^ //// 



+B{y) cos2(l)T 



(2h • pi^ h ■ 



P2t 



PIt ■ P2t 



Ml M2 



(5) 



where a is the fine structure constant, is the charge of a parton with flavor /, and 



A{y) 



(1 + cos^ e) 



B{y) = y(l -y)^= ^ sin^ i 



(6) 



The functions f(,h-^^ , are the distributions of unpolarized and transversely polarized partons with flavor / in 
unpolarized hadrons, respectively. Both are convoluted with their antiparton partners according to the definition 



fifl 



dpi^dp2T S {PIt +P2t - Qr) //(2:i,Plr)/l (2;2,P2t) + (1^2) 



(7) 



The cross section ((SJ can be integrated upon the angular distribution and the transverse momentum of the final 
lepton pair, 



da" 
dxidx2 



/7(xi)//(a:2) + (1^2) 



daf 



/7(xi)//(x2) + (1^2) 



(8) 



where, if neglecting the parton masses, s ~ X1X2S is the available cm energy for the elementary annihilation, and the 
elementary QPM cross section has been put into evidence. At this point, it is useful to recall that r = X1X2, and 
to introduce the new invariant x^ = xi — X2, which is the "longitudinal" momentum of the annihilating partons in 
their cm frame with respect to the maximum "longitudinal" momentum available. If da° is considered differential 
also in the leptonic invariant mass = M^, by transforming from the set {xi, X2} to the variables {r, x^}, and after 
integrating upon dr, we recover the well known result 



da" 



47ra^ 



dNPdx, 



Xl + X2 



E4 \f{ixi)f{ix2) + {1^2) 



(9) 



5 



about the scaling of the cross section in the cm energy s, which has been experimentally confirmed Wr\. Interestingly, 
when considered differential in -i/r and Xp, the cross section falls like ~ 1/s ; this remark will be reconsidered in the 
next Section because it influences the choice of the kinematical set-up. 



When one of the annihilating hadrons is transversely polarized, i.e. for the single-polarized Drell-Yan Hi H2 



IIX, 



the polarized part of the cross section becomes [34 



dAa^ 



E 



(Kldxidx2d(\^ 
-B{y) sin(3(/>-(/)5jJP 



h • pi 

( 



h p2, 



IT 



///: 



Mo 



Ml 

4h • Ply (h • p2y)^ - 2h • Pij 



P2j 



h Plr P2^ 



2Mi M| 



(10) 



where /-^j. is the Sivers function, describing the distribution of unpolarized partons with flavor / in a transversely 
polarized hadron, while h^rf! describes the distribution of transversely polarized partons with flavor / in transversely 
polarized hadrons. 

From Eq. H10() . a spin asymmetry weighted by sm{(j) — (/'S2) allows to isolate /f^ through the known distribution /i 
with a mechanism similar to the Sivers effect in DIS with lepton beams |i35] . Then, it could be possible to verify the 
interesting conjecture that the Sivers function should change sign when going from DIS to Drell-Yan because of 
the different behaviour of the tower operator that must be inserted in Eq. to link the space-time points and z 
in order to make the whole correlator color-gauge invariant. 

In Eq. Hl()(l . another interesting spin asymmetry is obtained by using the siii{(l) + (j)S2) weight, because the transver- 
sity hi can be extracted with a mechanism similar to the Collins effect However, the convolution involves also 
the unknown chiral-odd function h^, which luckily appears also in the leading- twist part of the unpolarized cross 
section (jSJ. A combined analysis of the cos 20 and sin(0 + (t>s2) moments of azimuthal asymmetries for unpolarized 
and single-polarized DrcU-Yan cross sections, respectively, should allow a complete determination of all the unknown 
distribution functions. The former asymmetry is particularly interesting because it could represent a natural expla- 
nation of the observed large azimuthal asymmetry of the cross section for the process 7r~ A —^ /i"*" fj,~ X , where A is 
either deuterium or tungsten, and the 7r~ beam energy ranges from 140, 194, 286 GeV [^SE^ 252 GeV 
leading order in as, the data can be parametrized as 



1 da" 
a° 'dn 



dxidx2dcij- 



J dxidx2dci^ T f( f( 

3 1/ 2 2 ^ i; 

; 1 + A cos 9 + u sin 6 cos (h -\ — sin 9 cos 2(h 

47rA + 3V ' ^2 



da° 



dndxidx2dq^ 



(11) 



with A 1 and /X < - 30 % |3a|39|. Both Leading Order (LO) and Next-to-Leading Order (NLO) perturbative 
QCD calculations give A ~ 1, but yit ~ ;/ ~ Other mechanisms, like higher twists or factorization breaking 

terms at NLO, are not able to explain the size of in a consistent picture |43 . l4al43| . On the contrary, from Eq. ||SJ) it 
is evident that a leading-twist (hence, large) azimuthal cos 2(j3 asymmetry in the unpolarized cross section is possible 
thanks to the distribution hj- [SJI- Since the latter is defined to describe the transverse polarization of quarks in a 
unpolarized hadron, it is intimately connected to the orbital motion of the parton inside the hadron; the product 
hj^ hi should bring a change of two units in the orbital angular momentum, such that the angular dependence involves 
twice the azimuthal angle (p. 

In the next Section, we will see how to numerically simulate these azimuthal asymmetries with a suitable Monte- 
Carlo. 



III. MONTE-CARLO SIMULATIONS 



In this Section we discuss a Monte-Carlo simulation for the specific processes pp fi'^ X and pp^ ^+ ^~ X. 
The primary goal of the simulation is to estimate the number of events required for such a Drell-Yan experiment in 
order to access unambiguous information on the distribution functions of interest, namely the transversity hi and 
hj^. Secondly, it is important also to clarify the role and the consequences of the unavoidable approximations and 
kinematical cuts. 
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The cross sections for the unpolarized and polarized processes, Eqs. (0 and (|10(l . respectively, are valid only in the 
framework of the QPM. In order to implement QCD corrections (and to test the resulting formulae), we will heavily 
refer to conventions and results of Ref. f4D], where a thorough analysis was performed on a large sample of muon 
pairs produced by a 252-GeV n~ beam interacting in a tungsten target. We will also use results from Ref. and 
Ref. 34] in adapting the relations of Ref. j40j to the case of an antiproton beam (with some care in matching the 
different notations) and of a polarized target, respectively. 

The Monte-Carlo events were generated from the following general expression for the cross section: 



da 



dndXp drdq^ 

or, equivalently. 



I 4 
^K-A{q^,x^,M)F{x^,T) c,(|q^|,a;^,r)5,(^,(/),0sj , (12) 



K^A{(i^,Xp-Xp,M){xp + Xp)F'{xp,Xp) ^ c',{\q^\, Xp, Xp) 3^(9, (j), (1)3^) , (13) 



dildxpdxpdq^ s 



=1 



where the above mentioned scaling in 1/s has been put into evidence. As the reader can argue from the formulae 
of previous Section, a (QPM-inspired) factorized structure for the cross section, but still differential in dq^ , implies 
that some simplifying assumption has been made on the dependence of the distribution functions, such that the 
convolution Q could be solved. The actual expression and its kinematical justification, as well as each other factor 
in the above equations, are separately discussed in the following. 



A. Kinematics and the transverse-momentum dependence 

For the handbag mechanism of Fig.^to be the dominant contribution, some kinematical contraints must be fulfilled. 
First of all, the cm energy s must be much bigger than the hadron masses involved. Secondly, in order to have a 
good control of the elementary annihilation process the invariant mass M of the muon pair must not overlap with the 
hadronic resonance spectrum. The most natural choice is to select M such that 4 GeV < M < 9 GeV: in this range, 
between the ip' and the first resonance of the T family, there should be no ambiguity in assuming that the elementary 
annihilation is followed by the creation of a virtual photon. 

We will consider both options where the antiproton beam hits a fixed proton target or it collides on a proton beam. 
In the first case, if the target is at rest, we have 

M'^ ^Ts^T [Pp + Ppf = T 2Mp{Mp + Ep)^T 2MpEp ; (14) 

near the limit t 1 (i.e., assuming that all the available cm energy flows into the final state), the upper limit for 
M would require a maximum antiproton beam energy of roughly 40 GeV. With the lower cutoff of M > 4 GeV, 
this corresponds to span the range 0.5 < a/t < 1, as discussed in the next Section and displayed in Fig. |3| Lower 
antiproton beam energies would reduce this scatter plot to an unrealistic small portion of phase space limited to very 
high values of both Xp and Xp, unless releasing the constraint on the lower cutoff of the dilepton mass. In fact, for a 
beam energy of Ep — 15 GeV it is also possible to cover a wide range of r with the constraint 1.5 GeV < M < 2.5 
GeV, which grants no overlap with the positions of the (/) and J/iJj resonances. Therefore, in the following we will use 
both energies for the antiproton beam. 

Recently, the discussions about the optimal setup configuration for spin asymmetry measurements at the HESR at 
GSI have focussed also on the collider option where two beams of protons and antiprotons collide at much higher s, 
particularly in the socalled asymmetric mode with Ep ^ Ep. Neglecting the hadron masses and in the same conditions 
as in Eq. l(T^ . we have 

^ {Pp + Pp f ^ 4EpEp , (15) 

because now Pp — —Pp. Hence, for the case of an antiproton beam energy of 15 GeV a proton beam of 3.3 GeV 
approximately gives s ~ 200 GeV^. By keeping the above cutoffs for M, a large portion of the phase space for Xp 
and Xp can be explored, particularly at lower values. Moreover, at larger energies the elementary process should 
be less affected by higher-order corrections like subleading twists; hence, the theoretical description should be more 
clean. Finally, as discussed below in Sec. IIII Dl the collider mode offers also the advantage that no dilution factors 
are introduced that affect the target polarization. For all these reasons, in the following we will present simulations 
for the single-polarized Drell-Yan process also at this kinematics. 
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The independent invariant fraction, 



Xp , can be redefined for practical purposes as [4C 



(16) 



where indicates the longitudinal component of the virtual photon momentum (in the following, we will also use 
the notation — jq^ |). The definition Hlt)|l puts in explicit evidence the content of x^ being the fraction of the total 
available longitudinal momentum in the collision cm frame. In the DIS regime, it recovers the previous definition in 
terms of light-cone momentum fractions. It is important to note that the minimum and maximum x^ values depend 
on the energy. For the case of Ep — 40 GeV and a fixed proton target, it turns out that —0.7 < a;^ < 0.7. 

A factorized transverse-momentum dependence in the differential cross section is achievable assuming a simple 
parametrization for the dependence of the distribution functions (typically, Gaussian- like) , which allows to directly 
calculate the convolution Q and get a factorized product of two functions depending separately upon and x^ , r 
(or, equivalently, Xp,Xp ; see, for example, the discussion of Sec. V in Ref. |3J|). Alternatively, since the input 
dependence and, consequently, the obtained q^, dependence are purely phenomenological, a parametrization for the 
latter can be assumed ab initio to fit the data. For example, the following form has been used in Sec. V of Ref. [iof. 
e.g. 



A(q^,X^,M) 



a 


It' 


a-1 


'b 


. b . 





(17) 



where a(xp,M), b{xp,M), are parametric polynomials given in App. A of the same reference. The normalization 
condition is 



/ 



dq^. A{q^ , Xf, , M) = 1 . 



(18) 



The function H17|) gives a good reproduction of the observed q^ spectrum for tt — p Drell-Yan events (see fig. 23 of 
Ref. 113] )■ However, in Ref. the same analysis has been repeated also for p — p events and the two distributions 
are similar for Q'^ > 1 GeV/c, while above 3 GeV/c big error bars prevent from any comparison. Therefore, we have 
adopted the same distribution ((T7jl also in our analysis. 

With the parametrization displayed in App. A of Ref. |4C|. a and b turn out to be smoothly dependent on x^ and 
Af, while the calculated (q^) is a smooth function only in M (after integrating on the whole range of x^^): it displays 
a strong decrease in x^. for x^. 0.7 after averaging on the whole range of M . The smooth dependence of (q^) on 
M reflects in its experimentally observed scale dependence ~ \/s, but if the cut-offs in the invariant mass are the 
same for two experiments at different energies (as it is the case between the measurement of Ref. ^40i] and the present 
simulation), the average squared transverse momentum results approximately the same. 

We have cut the distribution from below, because events at low q^ do not contribute to the spin asymmetry, but 
rather to the background, and they would represent a sort of artificial dilution factor. However, since this cutoff has 
drastic consequences on the final number of available events, a good compromise at these energies is represented by 
q^ > \ GeV/c, in accordance with several previous experimental observations. Such a constraint implies a reduction 
of the total number of available events by approximately 50%. 



B. QCD corrections and the partonic-momentum dependence 



The well established experimental observation that Drell-Yan pairs are distributed with {q^) > 1 GeV/c and 
depending on y/s, suggests that sizeable QCD corrections are needed on top of the QPM, because confinement alone 
induces much smaller quark intrinsic transverse momenta. The Feynman diagrams typically involve qq annihilations 
into gluons or quark-gluon scattering (see Fig. 3 of Ref. ^3)- They have been calculated introducing two main 
levels of approximation [sol l5ll | . The first one is the so-called Leading-Log Approximation (LLA) , where the leading 
logarithmic corrections to the Drell-Yan cross section can be resummed at any order in the strong coupling constant 
as, with the final net effect that the parton distribution functions get an additional scale dependence on M^. Following 
the prescriptions of Ref. 0] , the DGLAP evolution could be obtained by describing the functions with parameters 
explicitly depending on logM^ (see Apps. A,B and D in Ref. fo^' further details). However, since the M range 
explored in Ref. ^0] is practically the same assumed here, also the same parametrization described there in App. A 
is here retained; correspondingly, the function F' in Eq. (|13|) does no longer depend on the scale M^. We have 
explicitly verified that starting from the parametrization of Ref. '40| (and after a suitable transformation due to 
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FIG. 3: The scatter plot for 48000 events of Drell-Yan muon pairs produced by the colhsion of a 40 GeV antiproton beam on 
a proton target in the kinematic conditions discussed in the text. 

ecs I 

180^ -: ■ 1 

160 
140 
120 
100 

80 

60 
40 
20 

"o 20 40 60 80 100 120 140 160 180 

e Lab 



FIG. 4: The correlation between the polar angle of /i"*" in the laboratory, Siab, and in the CoUins-Soper frame, 9, for the collision 
of antiprotons with 40 GeV energy and fixed proton targets. 



different conventions), our Monte-Carlo gives a good approximation of the results of Ref. |43 with an antiproton 
beam of 125 GeV. 

The function F' in Eq. I|13() is given by 

2 

F'{xp, ^p) = ^Y. 4 // (^p) //(^p) + {P-^P)^ (19) 

namely it is the symmetric part of the unpolarized cross section that has been factorized out for convenience (see 
Eq. ©). With fi {x) parametrized as in App. A of Ref. for the various flavors f = u,d, s, we show in Fig. O the 
scatter plot in Xp and Xp of some 48000 events sorted by the cross section of Eq. (|13ll including Eq. H19|l and assuming 
the kinematics corresponding to Ep = 40 GeV with the cutoffs > 1 GeV/c, 4 < M < 9 GeV and -0.7 <Xp< 0.7. 

The line bisecting the plot at 45 deg. corresponds to x^ = 0; parallel lines above and below it indicate a;^ > 
and Xp < 0, respectively. Hyperboles with XpXp — const, select various values of r: moving to the upper right 
corner of the figure corresponds to the limit t — > 1, where events are rare; on the contrary, approaching the lower 
cutoff -^/r ~ 0.5 the distribution becomes much more dense. Nevertheless, the projection along each axis displays a 
significant number of events also at the boundary values of Xp,Xp. They correspond to large positive and negative 
values of Xp , which in turn indicate, in the CoUins-Soper frame, small and large values of the muon pair polar angle 
9, respectively. For a fixed target and in the laboratory frame, since we have antipartons in the beam and partons 
in the target, large and positive values of x^, correspond to "forward" events where, by convention, the polar angle 
^lab of fi'^ is small. Viceversa, large negative values of Xp correspond to "backward" events with large ^lab. This is 
at variance with the experiments of Refs. where events with only large positive x^ were collected because of 

the higher beam energy and of the stronger absorption of the unpolarized target. In Fig. ^ the correlation between 
9 and Oiah is shown in order to estimate the possible geometry of the apparatus for detecting the muons. The points 
are spread because the relation between the two frames depends on Xp,Xp, and q^. 

The second level of approximation in the QCD corrections to the QPM result is referred to as Next-to-Leading-Log 
Approximation (NLLA). It amounts to include in the calculation all processes at first order in as involving a quark. 
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an antiquark and a gluon |50l l51|. The effect is large, because the cross section has approximately twice the size of 
the corresponding one in the QPM. This correction turns out to be roughly independent of a:^ and (except for 
the kinematical upper limits) and it is usually indicated as the K factor. As explained in Ref. (49], K depends on the 
choice of the parametrization of the distribution functions through their normalization, but it also grows like ^/t |4Clj| . 
However, in the following we will consider azimuthal asymmetries, which are defined as ratios of cross sections. Hence, 
it is legitimate to expect that such corrections in the numerator and in the denominator approximately compensate 
each other. Indeed, for fully polarized Drell-Yan reactions at high energy the spin asymmetry is weakly affected by 
QCD corrections j23|. In any case, for the considered antiproton energy Ep = 40 GeV at ~ 0.5, where most of 
the events are concentrated, the best compromise is to assume the constant value K = 2.5. 



C. The angular distribution 

If the orientation of the x axis of the CoUins-Soper frame (pointing along ) with respect to the laboratory frame 
is indicated with then it turns out that 0' — —4>Sp because, in practice, the latter angle gives the position of the 
X axis of the laboratory frame (pointing along Sp^) in the CoUins-Soper frame. The whole solid angle (^,0) of the 
final muon pair in the Collins-Soper frame is randomly distributed in each variable. The explicit form for sorting the 
angular distribution in the Monte-Carlo is 



4 

E 

i=l 



c[{q^,Xp,Xp) Si{9,(l), (f)Sp) 



1 



+ cos^ 9 

v(xp, Xp, Qj.) 



sin 6 cos2(/) 



+ \Sp^\c4{q^,xp,Xp) Si{9,(l),(j)Sj,) , 



(20) 



which corresponds to the parametrization for A = = 0, as suggested by experimental observations, and 
includes also a polarized term as suggested by Eq. l(TT)|l . 

If quarks were massless, the virtual photon would be only transversely polarized and the angular dependence would 
be described by the functions ci = = 1 and C2 = 1, 5*2 = cos^ 9. Violations of such azimuthal symmetry induced 
by the function C3 = are due to the longitudinal polarization of the virtual photon and to the fact that quarks have 
an intrinsic transverse momentum distribution. QCD corrections influence such function, which in principle depends 
also on Xp,Xp, and (see App. A of Ref. ^Ul)- practice, we follow Ref. [s^ (where v — 2k and k is given in 
Sec. V) and assume the simple form 



v{q^) = 3.5 



4M2q2 



(4M2 



,2 ^2 



(21) 



where = 2.3 GeV. 



The last term in Eq. (|20|l corresponds to the polarized part of the cross section, which is fully described in Eq. I|10|l . 
Since we are interested in the convolution of hi and hj^, we assume that the angular dependence of the azimuthal 
spin asymmetry is given by 



^4 



J 



sm 



J- 



(22) 



Recalling that in Eq. H13|) the azimuthally symmetric unpolarized part F'{xp,Xp) of the cross section has been 
factorized out, the corresponding coefficient C4 in Eq. (|20|) becomes 











[^■P^- Mp \ 









(23) 



Following Ref. [34j . we introduce some model assumptions on the behaviour of the unknown /ij'-(x,Py). Its 
dependence is of the Gaussian type, while the x dependence is given directly by fi{x): 



TT p2 



Ml 



fiix) , 



(24) 
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where = 1 GeV ^. Then, the C4 coefficient turns out to be (see Sec. VI of Ref. [s^l for all details): 

cJa XX)- Efe}f{{xp)h{ix,) 

4iW^ + 9t L,/ 6/ /i {Xp) f( {Xp) 

For sake of simplicity, we assume that we can approximate the contribution of each flavor to the parton distributions 
by a corresponding average function, such that we finally get 

^ 2g^^c {fixp)){hi{xp)) _ 2q^Mc {hijxp)) 

c,iq,,xp,xp)^ 4M2 +,2 (/>,))(/,(:,,)) - 4Af^+gJ (A(xp)) • ^''^ 

In order to disentangle the various azimuthal parts in the cross section, the most straightforward way is to inte- 
grate the corresponding azimuthal asymmetry with a proper weight depending on (f). However, if the sin(3(^ — (t>S2) 
contribution to Eq. pOf) is assumed small, a simpler procedure is available. In fact, for the remaining leading-twist 
contributions 

P;^(0) = cos2(/) , 
P2(0) = sin(0-05j , 

Psicb) = sin(0 + (/.5j , (27) 

it is possible to show that for each bin Xp the asymmetry between the subsets of positive and negative values of Pi{(f>) 
(spanning the whole (j) range) does not depend on Pj for j 7^ i. In the following Section, results will be presented and 
commented for azimuthal asymmetries constructed in this way for Pi{4>) and P^{(j)). Two choices with opposite features 
will be selected for the ratio {hi{xp)) / {fi{xp)) , namely the ascending function 2xp and the descending one 2(1 — Xp). 
The goal is to determine the minimum number of events (compatible with the kinematical setup and cutoffs) required 
to produce azimuthal asymmetries that can be clearly distinguished like the corresponding originating distributions. 
In fact, this would be equivalent to state that some analytic information on hi(x) could be extracted from such spin 
asymmetry measurements. 



D. Polarized target and dilution factor 

When the transversely polarized proton is considered as a fixed target, in reality we assume a NH3 target with 85% 
transverse polarization for the H nucleus. In the ideal conditions of scattering from a black sphere, the N hadronic 
cross section is exactly ten times bigger than the one for H, i.e. 390 mb ^5^. This means that, from the point of view 
of hadronic collisions, the N nucleus is equivalent to 5 protons and 5 neutrons with no difference between protons and 
neutrons. In reality, this picture can be valid only at very small Xp and Xp, where the valence quarks do not dominate. 
However, as it is clear from Fig.Othis domain is not accessible to our simulation. Rather, the events are concentrated 
at Xi > 0.1, i = p^p, where the valence distributions are dominant and the annihilation occurs between a quark and 
an antiquark of the same flavor. Hence, the antiproton of the beam "prefers" to hit protons rather than neutrons. 
A convenient estimate would be to assume that in the NH3 target the N nucleus behaves like if it was made of 5 
protons and 4 neutrons, such that out of 12 events, 9 come from collisions inside N and 3 inside H3. Therefore, in the 
scatter plot of Fig. O with 48000 events, only 12000 are related to with 85% of transverse polarization, while 36000 
are produced on the unpolarized N nucleus. The Monte-Carlo simulation takes into account such a combination of 
dilution factors. 

When the pp^ jX^ X process is considered as the collision of two beams, a systematic dilution factor 0.85 is 
applied to the proton polarization. The statistics of the simulation greatly benefits from this weaker dilution and, in 
fact, a much smaller sample of events is needed to get reasonable results, as it will become evident in the next Section. 



IV. RESULTS 

In this Section, we present results of a Monte-Carlo simulation for the pp^ — > ji^ /i^ X and pp /i+ /i^ X processes 
in order to make realistic estimates of the minimum number of events required to extract as detailed information as 
possible on the chiral-odd distributions hi and h^. First, we consider an antiproton beam of 40 GeV energy hitting 
a fixed proton target and producing muon pairs with invariant mass in the range 4 GeV < Af < 9 GeV and total 
transverse momentum > 1 GeV/c. Then, we will select the lower antiproton beam energy of 15 GeV with the 
dilepton invariant mass in the range 1.5 GeV < M < 2.5 GeV and we will explore both possibilities of a fixed proton 



11 



TABLE L Total absorption cross sections per nucleon for antiprotons at the indicated energies and for various invariant masses 
of the Drell-Yan pair (see text for a discussion of the kinematics and the cutoffs). 



Ep (GeV) mode event sample M (GeV) atot (nb/nucleon) 

40 fixed target 40000 4^9 0.02 

40 fixed target 40000 1.5^2.5 1.6 

15 fixed target 40000 4 4 x 10"** 

15 fixed target 40000 1.5-^2.5 0.8 

15 colhder 8000 4^9 0.1 

15 colhder 8000 1.5-^2.5 2.4 



target and of a collision with a proton beam of 3.3 GeV energy, such that the cm available squared energy is s 200 
GeV^ (see Sec. lIII Al . For fixed targets at such low energy, the QPM picture leading to Eq. ((T^ . or ((T^ . is not well 
justified because of the unavoidable corrections due to higher-twist contributions. But we stress that the primary goal 
of this simulation is to explore the feasibility of the extraction of the hi and hj^ distributions, not to make a precision 
calculation. Hence, we will retain the same approach in all kinematical cases above described. 

The transversely polarized proton target is obtained from a NH^ molecule where each H nucleus is 85% transversely 
polarized and the number of "polarized" collisions with each H nucleus is 25% of the total number of collisions. The 
events are sorted according to the cross section of Eq. (|13|l (reduced by the above dilution factor) and supplemented 
by Eqs. H17|l and H19l) - (|26|) . In the following, we will study the azimuthal asymmetries generated by the cos 2^ and 
sin((/) + (l>Sp) dependences in the cross section (see Eqs. ©, (flUll . (1^ and (E^J). For each Xp bin, the asymmetry 
is constructed by dividing the sample in two groups, one for positive values in the dependence (U) and another 
one for negative values (D), and taking the ratio {U — D)/{U + D); the azimuthal cj) dependence is considered in the 
CoUins-Soper frame. 

We have initially taken samples of 100000 and 20000 events for fixed target and collider modes, respectively. The 
above cuts in M and are implemented also by a restriction on the 9 angular dependence to the range 60 deg. 
< 9 < 120 deg., where events are mostly concentrated; outside these limits, in fact, the azimuthal asymmetries 
are too small. Analogously, as a general rule Xp bins with less than 100 events are considered empty. All these 
cuts reduce the starting samples approximately by a factor 2.5; conventionally, we will indicate the surviving 40000 
and 8000 events, respectively, as "good" events. Statistical errors for the above mentioned azimuthal asymmetry 
{U — D)/{U + D) are obtained by making 20 independent repetitions of the simulation for each considered case, 
and then calculating for each Xp bin the average asymmetry value and the variance. We checked that 20 repetitions 
are a reasonable threshold to have stable numbers, since the results do not change significantly when increasing the 
repetitions from 6 to 20. 

From our Monte-Carlo results, we deduce a total cross section atot — 0.02 nb per single nucleon for the absorption 
of antiprotons in a 40 GeV beam hitting a polarized NH3 target in the above discussed conditions and producing 
final lepton pairs with invariant mass M > 4 GeV. This cross section is 2 x 10^ times smaller than the total nucleon 
absorption cross section per single nucleon at the considered energy and under the conditions for a complete absorption, 
since the latter approximately amounts to 40 mb. In other words, when the flux of actually absorbed antiprotons is 
deduced from the experimental setup, the above cross section states that one out of 2 x 10^ such antiprotons produces 
a "good" event for the selected kinematics. In Tab.Q] we list the total cross sections for all the explored kinematics and 
for various ranges of invariant masses. If C is the luminosity in units (cm~^s"^), the product <Jtot^ gives the number 
of "good" events per nucleon and per second. For example, we deduce that in the collider mode for 1.5 < M < 2.5 
GeV we would approximately get 60000 "good" events/month with C = 10'^^ (cm"^s~^). 

A. Polarized Drell-Yan 

In Fig. [SJ the sample of 40000 "good" events for the pp^ fi'^ fi~ X process is displayed for the antiproton 
beam energy Ep — 40 GeV and for the case of a fixed proton target. The left panel corresponds to the choice 
{hi{xp)) / {fi{xp)) = 2xp, while the right one to the function 2(1 — Xp), according to Eq. Results are reported in 

Xp bins over the entire range, but the bins at the boundaries, corresponding to Xp < 0.2 and Xp > 0.8, are scarcely or 
not at all populated, according to the phase space distribution of Fig. 13 For each bin, two groups of events are stored: 
one corresponding to positive values of sin(0 + (ps ) in Eqs. (|20|l and H22|) (represented by the darker histogram), the 
other one corresponding to negative values (represented by the superimposed lighter histogram). 

In Fig. the resulting asymmetry {U — D)/{U + D), between cross sections with positive (U) and negative 
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FIG. 5: The sample of 40000 events for the pp^ ^ fi'^ ii~ X process on a fixed proton target with an antiproton beam energy 
Ep — 40 GeV and with a muon pair of invariant mass 4 < M < 9 GeV and transverse momentum q^. > 1 GeV/c (for 
further details on the cutoffs, see text), a) Left panel for the choice {hi{xp)) / {fi{xp)) = 2xp (brackets mean that each flavor 
contribution in the numerator is replaced by a common average term, similarly in the denominator; for further details, see 
text), b) Right panel for {hi{xp)) / {fi{xp)) = 2(1 — Xp). For each bin, the darker histogram corresponds to positive values of 
sin(<^ + cpSp) in Eqs. 12011 and 1221 . the superimposed lighter one to negative values. 
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FIG. 6: Azimuthal asymmetry {U — D)/{U + D) between cross sections in the previous figure corresponding to the 
darker histograms (U) and superimposed lighter histograms (D), as bins in Xp. Downward triangles for the case when 
{hi{xp)) / {fi{xp)) = 2xp, squares when it equals 2(1 — Xp). Continuous lines are drawn to guide the eye. Error bars due 
to statistical errors only, obtained by 20 independent repetitions of the simulation (see text for further details). 

(D) values of sm{(j> + (j)Sp), is plotted again in Xp bins for both choices: downward triangles for the case when 
{hi{xp)) / {fi{xp)) = 2xp, squares when it equals 2(1 — Xp). The error bars represent statistical errors only. 

Here and in the following cases, both functional forms for {hi{xp)) / {fi{xp)) , once integrated in Xp over the entire 
range, give the same overall size 1. Extrapolating information from other sources (lattice calculations, other exper- 
iments like semi-inclusive DIS, theoretical properties like the Soffer bound, etc..) suggests that this size may be a 
reasonable expectation js^l- From Fig.|H| we first conclude that 40000 "good" events of antiprotons with energy 40 
GeV hitting fixed protons and producing muon pairs with 4 < Af < 9 GeV, are enough to produce a statistically 
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FIG. 7: Same as in Fig. 13 but witli the antiproton beam energy Ep = 15 GeV and with a muon pair of invariant mass 
1.5 <M< 2.5 GeV. 
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FIG. 8: Azimuthal asymmetry coUected in Xp bins in the same kinematics as in Fig.Qand with the same notations as in Fig. El 
Continuous hues are drawn to guide the eye. Error bars from statistical errors only. 

nonvanishing azimuthal asymmetry of about 5% on the average. But it is difhcult to disentangle the two different 
functional forms of the ratio except for very high Xp > 0.6, where the QPM picture becomes questionable. Equiva- 
lently, it is difficult to extract unambiguous information on the transversity distribution hi{x). Moreover, from Tab.Q 
we deduce that a luminosity C — 10^^ cm~^s~^ would give a rate of 500 "good" events/month. The same conclusion 
holds if the absolute size of the ratio is expected in the range 0.5-;- 1. Only a smaller number by an order of magnitude 
would imply a much larger sample of events in order to reach the same sensitivity. 

The situation is somewhat better at the lower antiproton beam energy Ep = 15 GeV, because in the scatter plot 
(analogous of Fig. |3Jl the lowest cutoff in r = XpXp = M'^/s is smaller since now 1.5 < M < 2.5 GeV. The highest 
density of events (and the best statistics) is now reached at lower values of Xp. In fact, while in Fig. [T] we show the 
usual sample of 40000 "good" events in the same conditions and notations as in Fig. |Slbut for Ep = 15 GeV and 
1.5 < AI < 2.5 GeV, in Fig.|Slthe resulting azimuthal asymmetry (again, with the same notations as in Fig.EJ clearly 
shows how the ascending (descending) behaviour of {hi{xp)) / {fi{xp)) = 2xp [2(1 — Xp)] reflects in the corresponding 
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asymmetry making the two results statistically distinguishable for Xp < 0.4 (and, approximately, also for very high 
values Xp > 0.7). 
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FIG. 9: The sample of 8000 events for the pp' —> /j. ^~ X process where an antiproton beam with energy Ep — 15 GeV 
collides on a proton beam with energy Ep = 3.3 GeV producing muon pairs with invariant mass 4 < M < 9 GeV and 
transverse momentum > 1 GeV/c (for further details on the cutoffs, see text). The content of each panel and the notations 
are as in Fig. |3 

The lowest range in the muon pair invariant mass offers also a faster rate in the collection of "good" events. From 
Tab. m we deduce that with the above mentioned luminosity C — 10'^^ cm^^s^^ it should be possible to collect a 
sample 40 times bigger than in the previous case, namely 20000 "good" events/month; or, equivalently, to reach the 
same counts (500 "good" events/month) with C = 2.5 x 10^^ cm~^s~^. 

From previous examples, it seems evident that the smaller the lower cutoff in r, the lower the range in Xp where 
the highest density of events gives the best available statistics to extract as many details as possible about hi. Since 
T = M'^/s and we want no overlap with M values corresponding to hadronic resonances whose elementary mechanism 
for the production of the final muon pair is not well under control, it might be convenient to considerably increase s. 
The most effective way is to consider the proton not as a fixed target but as a beam colliding against the antiproton 
one. Keeping for the latter the energy Ep — 15 GeV, a proton beam with Ep = 3.3 GeV produces an available 
cm squared energy s — 200 GeV^ (see Sec. IIII Al for details). In Fig. |51 a sample of 8000 "good" events for the 
pp^ — *■ /Lt^ /x~ X process leading to muon pairs with invariant mass 4 < M < 9 GeV is plotted in Xp bins, all the 
notations and the other cutoffs being the same as before. Evidently, events are concentrated at lower Xp bins; but 
also the asymmetry is, as it can be deduced from the different color codes of the histograms. In fact, in Fig. ^| 
the corresponding azimuthal asymmetry ([/ — D)/{U + D) is plotted again with the same notations as in Fig. El 
the Xp bins being limited to 0.6 because the histograms corresponding to higher Xp values are empty. First of all, 
the absolute size of such asymmetry is much bigger (note the different scale with respect to Figs. El and |H1); on the 
average, it is approximately 15%. But, most important, the best statistics is achieved for Xp < 0.4, where the small 
error bars allow for a very clean distinction between the results corresponding to the opposite trends for the function 
{hi{xp)) I {fi{xp)) . And all this can be obtained with a reduced sample of just 8000 "good" events. 

In conclusion, although the option of a fixed proton target permits a sufficiently clean extraction of the analytic 
properties of the function hi{x) for x < 0.4 provided that the invariant mass of the final muon pair is enough low, 
the collider mode seems more promising. The higher available cm energy grants a cleaner theoretical description of 
the elementary mechanism, gives bigger azimuthal asymmetries, and allows for a much better statistics and a clean 
information in the same range x < 0.4. As for the event rate, from Tab. jHwe deduce that the luminosity C = 10'^^ 
cm~^ s~^ would allow to collect 2500 "good" events/month. 
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FIG. 10: Azimuthal asymmetry collected in Xp bins in the same kinematics as in Fig. El and with the same notations as in 
Fig-ini Continuous lines are drawn to guide the eye. Error bars from statistical errors only. 
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FIG. 11: The sample of 40000 events for the unpolarized pp — > /i^ /i~ X process on a fixed NHz target with an antiproton 
beam energy Ep = 15 GeV and with a muon pair of invariant mass 1.5 < M < 2.5 GeV (for further details on the cutoffs, 
see text), a) Left panel for 1 < Pj, < 2 GeV/c. b) Right panel for 2 < p^ < 3 GeV/c. For each bin, the darker histogram 
corresponds to positive values of cos2<^ in Eq. 11201 . the superimposed lighter one to negative values. 



B. Unpolarized Drell-Yan 

In Fig. ^2 the sample of 40000 "good" events for the unpolarized pp ^ fi^ fi~ X process is displayed for the case of 
a fixed NH^ target and for the antiproton beam energy Ep = 15 GeV, i.e. in the same initial conditions proposed at 
the HESR at GSI by the PANDA collaboration [H^. The muon pair invariant mass is constrained by 1.5 < M < 2.5 
GeV and 60 deg. < 9 < 120 deg. The left panel corresponds to muon pairs with transverse momentum 1 < < 2 
GeV/c, while the right one to the choice 2 < p^ < 3 GeV/c. The goal is to test the simple distribution assumed 
for the coefficient i/ in Eq. (|21|l which leads to the corresponding factorized contribution in Eq. H2U|) . For each bin, 
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two groups of events are stored corresponding to positive (darker histograms, U) and negative (superimposed lighter 
histograms, D) values of cos 2^ in Eq. (|20|) . Again, as in previous plots results are reported in Xp bins over the entire 
range, but the bins at the boundaries, corresponding to Xp < 0.2 and Xp > 0.8, are scarcely or not at all populated 
and they will be discarded in the plot of the corresponding azimuthal asymmetry {U ~ D)/{U + D). 
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FIG. 12: Azimuthal asymmetry {U — D)/{U + D) between cross sections in the previous figure corresponding to the darker 
histograms (f/) and superimposed lighter histograms (D), as bins in Xp. Downward triangles for muon pairs with transverse 
momentum 1 < p^, < 2 GeV/c, squares for 2 < < 3 GeV/c. Continuous lines are drawn to guide the eye. Error bars due to 
statistical errors only. 

In Fig. 1121 the resulting asymmetry {U — D)/{U + D) between cross sections with positive iU) and negative {D) 
values of cos 20 in Eq. (|20|l , is sorted in Xp . This asymmetry, already observed at higher energy [4fij , is responsible 
for the violation of the so-called Lam- Tung sum rule and, according to Eq. (O, it could be related to the distribution 
function hj^ , namely it could be interpreted as a distortion of the parton momentum distribution due to its transverse 
polarization. The predictions of Fig. I12| should be reliable, because the function v{q^) has been tested against the 
most recent available measurements 40]. In fact, despite its simple functional dependence viq^,) should not depend on 
the energy; therefore, even with an antiproton beam of 15 GeV the expected size of the asymmetry should be the one 
displayed in Fig. 1121 Moreover, in the unpolarized Drell-Yan the magnitude of the asymmetry should also not depend 
on the target, because there is no dilution factor as it is the case for the corresponding polarized process. Nevertheless, 
in Fig. ll2l two different sets of asymmetries are plotted corresponding to different ranges for the transverse momentum 
of the final muon pair: downward triangles for 1 < < 2 GeV/c and squares for 2 < < 3 GeV/c. The good 
statistics and the small error bars should allow to reliably test the q^^ dependence assumed for the coefficient v in 
Eq. laj. 

V. CONCLUSIONS 

In this paper, we have concentrated on the investigation of the spin structure of the proton using the single-polarized 
Drell-Yan process pp^ At leading twist, the polarized part of the cross section contains three contributions 

with different dependences on the azimuthal angle (f> identifying the direction of the final muon pair [s^ . We have 
focussed on the one involving the convolution of the transversity distribution, the missing piece necessary to complete 
the knowledge of the nucleon spin structure at leading twist [a, 0] , and the distribution hj^ , which describes how in 
an unpolarized hadron the momentum distribution of a parton is distorted by its transverse polarization. Extraction 
of the latter is of great importance as well, because hj~ contains information on the orbital angular momentum of 
partons and it appears also at leading twist in the unpolarized part of the cross section,^ weighted by cos 2(j). As such, 
it is believed to be responsible for the well known violation of the Lam- Tung sum rule 39, "4^, an anomalous big 
azimuthal asymmetry of the unpolarized Drell-Yan cross section, that neither NLO QCD calculations ^'^^ higher 
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twists or factorization-breaking terms in NLO QCD Q| are able to justify. Measuring at the same time the 

azimuthal asymmetries of the unpolarized and single-polarized Drell-Yan cross section allows to determine both the 
unknown functions, i.e. the transversity and hj^. 

Using unpolarized antiproton beams offers the advantage of involving (anti)parton valence distributions that should 
not be suppressed as in the standard fully polarized Drell-Yan |27i, ,28] , but avoiding at the same time the difficulty 
of polarizing antiprotons. At GSI a new facility is under development, the High Energy Storage Ring (HESR), where 
an antiproton source will be built. The transversely polarized proton target could be obtained from a NH3 molecule 
where each H nucleus, for example, is 85% transversely polarized; a further 25% dilution factor needs to be applied 
to filter out spurious unpolarized events from the N nucleus. We have simulated the above single-polarized Drell-Yan 
process using an antiproton beam of 40 GeV and selecting muon pairs with invariant mass M in the range 4 9 
GeV and with a total transverse momentum above 1 GeV/c. In this kinematics, the quark-antiquark pair annihilates 
into a virtual photon and the theoretical description of the elementary process is well under control. We have also 
made simulations at the lower beam energy of 15 GeV but with invariant masses in the range 1.5 2.5 GeV, in order 
to keep a reasonably large phase space and at the same time to avoid the overlap with hadronic resonances like the 
at 1 GeV and the J/ip at 3.1 GeV. For such low M values, the parton model interpretation of the elementary 
annihilation is questionable because of the contributions of higher twists. But from the analysis of the scatter plot 
events seem to concentrate at the lowest possible values of = M'^/s (where the fractions of 

hadron momenta carried by the annihilating quark-antiquark pair and s is the available initial squared energy in the 
center-of-mass system); hence, the simulation is still useful to study the gross features of the cross section and to 
explore the possibility of extracting the transversity distribution at this lower energy. For both beam energies, the 
polar-angle distribution of the muon pairs is restricted to the range 60 120 deg., where most of the Monte-Carlo 
events are concentrated. An alternative way to reach the high-statistics portion of the phase space at low r is to 
increase s. Therefore, we have considered also the option where the antiproton beam with energy Ep = 15 GeV 
collides on a proton beam with Ep = 3.3 GeV such that s — 200 GeV^. In this last case, 4 < M < 9 GeV, all the 
other applicable cutoffs being unchanged. 

The primary goal of the simulation is to estimate the number of events required in order to access unambiguous 
information on the distribution functions of interest. Azimuthal asymmetries have been built, for each bin in the 
partonic momentum fraction, by separating the cross sections with positive values in the (p dependence from the 
negative ones and taking the ratio between the difference and the sum of such values. The cos 2(/) asymmetry in the 
unpolarized cross section has already been measured 40] at higher energies. However, under the present kinematical 
conditions the magnitude of the asymmetry should be independent from the energy of the experiment and it can 
be used as a cross-check of the Monte-Carlo simulation. It turns out that with an initial sample of 40000 events 
(surviving the filter of the various discussed cutoffs) it is possible to study in detail the function v responsible for such 
an asymmetry, but in the restricted range 0.2 ~ 0.8 for the parton momentum fraction, because the boundary bins 
are not enough populated. 

In the case where antiprotons with energy 40 GeV hit transversely polarized fixed proton targets and produce muon 
pairs with 4 < M < 9 GeV, the 40000 events are enough to produce a statistically nonvanishing azimuthal asymmetry 
of about 5% on the average, but it seems difficult to extract more detailed information on the analytic structure 
of hi(x) except for very high x > 0.6, where the adopted QPM picture becomes questionable. The lower invariant 
mass range 1.5 < M < 2.5 GeV allows to reach lower values of r in the phase space, where most of the events are 
concentrated, even for lower antiproton beam energies such as 15 GeV. Hence, a better statistics is reached for lower 
Xp values and it seems possible both to extract the functional dependence of hi{x) for x < 0.4, and to collect events 
at a higher rate (several times bigger than the previous case at the same luminosity). However, the best option seems 
the one where antiproton and proton beams collide reaching high values of s. In the explored case s ~ 200 GeV^, 
the theoretical description of the elementary mechanism is clean and is not affected by complications like, e.g., higher 
twists; moreover, the absolute size of the azimuthal asymmetry turns out bigger and the much better statistics should 
allow for an unambiguous extraction of hi{x) for x < 0.4 with a reduced sample of just 8000 events. 

In conclusion, we hope that the present work will help in the feasibility studies about the physics program of 
hadronic collisions with antiproton beams at the HESR at GSI. 
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